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Abstract 

We give a Hamiltonian based interpretation of microscopic Fermi liquid 
theory within a renormalization group framework. We identify the fixed point 
Hamiltonian of Fermi liquid theory, with the leading order corrections, and 
show that this Hamiltonian in mean field theory gives the Landau phenomeno- 
logical theory. A renormalized perturbation theory is developed for calcula- 
tions beyond the Fermi liquid regime. We also briefly discuss the breakdown 
of Fermi liquid theory as it occurs in the Luttinger model, and the infinite 
dimensional Hubbard model at the Mott transition. 
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1 Introduction 



The recent theoretical attempts to understand the anomalous behaviour of the 
high temperature superconductors have raised some fundamental questions as to the 
form of the low energy excitations in strongly correlated low dimensional Fermi sys- 
tems. Before the discovery of the high T c materials the superconductivities in most 
superconducting metals and compounds have been explained within the BCS theory 
as instabilities within a Fermi liquid induced by an attractive retarded interaction 
due to coupling with the phonons. Possible exceptions are the heavy fermion super- 
conductors, such as UPt 3} where purely electronic mechanisms have been proposed 
for the effective interelectron attraction. Even in these cases the superconducting 
instability is believed to be one within a Fermi liquid. However, in their normal 
state high temperature superconductors are not good metals and their behaviour 
appears to differ from that of a conventional Fermi liquid (for a review of the ex- 
periments on these systems see [f]). This has led to conjectures that the normal 
state has no well defined quasi-particles at the Fermi level so that it can not be 
a Fermi liquid; there have been conjectures that the behaviour is better described 
by some form of 'marginal Fermi liquid' [2] or 'Luttinger liquid' [3]. More extreme 
breakdowns of Fermi liquid theory have also been proposed where no Fermi surface 
remains in the usual sense because the imaginary part of the self-energy is always 
finite [4]. These are still a controversial issues. However, these theories all propose 
that the superconductivity in the high T c materials is not an instability in a Fermi 
liquid and that we are dealing with a novel situation. As the characteristic feature 
of these materials is the CuO^ plane, the basic question is whether the Fermi liquid 
breaks down in these two dimensional systems due to strong correlations induced 
between the d electrons at the Cu sites. In one dimension it is well established that 
no matter how weak the inter-electronic interaction the Fermi liquid theory breaks 
down and the low energy excitations for short range repulsive interactions are well 
described by a suitably parametrized Luttinger liquid. Anderson's theory [3] for the 
high T c materials is based on the assumption that a similar Luttinger liquid state 
occurs in strongly correlated two dimensional models. If this is so, then could such 
a breakdown occur for higher dimensional systems with strong correlation? Is there 
a critical interaction strength for the breakdown of Fermi liquid theory for the two 
dimensional systems? 

These are the sort of questions have led to a re-examination of Fermi liquid 
theory. It is nearly forty years since Landau proposed his phenomenological Fermi 
liquid theory [5] , and it was very soon after that the approach was verified for mi- 
croscopic models within the framework of many-body perturbation theory [6] [8] [7]. 
Since that time new techniques for tackling problems in condensed matter physics 
have been devised, such as the renormalization group approach, and others, such 
as 'bosonization', have been developed further. These approaches can bring a fresh 
perspective to the subject. Whether any of these techniques can provide answers to 
the question of the behaviour of two dimensional Fermi systems cannot be answered 
at this stage. What is clear, however, is that they can give new insights which are 
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of value in themselves, and provide a new context for addressing these problems. 

The bosonization approach in one dimension gives the essential physics for many 
interacting fermion models in a very simple way. This was shown in the pioneering 
work of Tomonaga [9], Luttinger [10], Mattis and Lieb [11], then further developed 
by Luther [12], Emery [13], Haldane [14] amongst others. Haldane recognized uni- 
versal features in these results and showed that low energy behaviour of a large 
class of one dimensional interacting fermion systems can be derived from a suitably 
generalized Luttinger model. He coined the term 'Luttinger liquid' to describe these 
in analogy with concept of a Fermi liquid for higher dimensional Fermi systems. 
More recent bosonization work [15] [16] has been concerned with showing how this 
approach can be generalized to higher dimensions, and showing how the Fermi liquid 
theory can emerge from bosons when they are constrained to the region of the Fermi 
surface. In this article, however, we will focus our attention on the renormalization 
group approach, both as developed by Wilson [17] in the 70s to tackle problems of 
critical phenomena but also as originally developed in field theory, QED etc, as a re- 
organization of perturbation theory. Recently the Wilson renormalization approach 
to these fermion models has been the subject of a review by Shankar [18]. In this 
article we cover similar ground but with a different emphasis and different examples 
so that there is in detail little overlap between the two surveys. The basic message, 
however, I think is the same, that the renormalization group gives us new insights 
into Fermi liquid theory, its stability and the circumstances which may cause it to 
breakdown. We develop a renormalized perturbation theory which can be directly 
related to the Wilson type of calculations, and also to the microscopic derivation of 
Fermi liquid theory. 

Before embarking on the renormalization group approach we briefly look at the 
main features of Fermi liquid theory as originally introduced by Landau. The 
Landau phenomenological theory [5] is based on the assumption that the single 
(quasi)particle excitations at very low temperatures of an interacting Fermi system 
are in one-to-one correspondence with those of the non-interacting system. In terms 
of the one electron states a total energy functional E to t is constructed of the form, 

E t ot = E gs + Y^ ^a}j Ua ^ + \ f°,°'( a -> a')Sn a}(T Sn a >y + ... (1) 

where E gs is the ground state energy, 6n a ^ is the deviation in occupation number of 
the single particle state |a),cr with an excitation energy from its ground state 

value, and is the leading term due to the quasiparticle interactions. A free en- 
ergy functional F is constructed from (1), retaining only the first two terms, together 
with the Fermi-Dirac form for the entropy of the quasiparticles. Minimization of 
F with respect to 8n a leads to asymptotically exact results for the thermodynamic 
behaviour as T, H — > for systems in a normal paramagnetic ground state. The 
reason why the higher order terms in (1) give negligible effects in this limit is be- 
cause the expectation value of Sn a , (Sn a ) — > as T — > and H — > 0. The effective 
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quasiparticle energy e„ j(7 in the presence of other excitations is given by 

£a,a = 4°t + fa,a'(a,a')(Sn a fy). (2) 

a',a' 

From this result asymptotically exact results for the specific, susceptibility and other 
low temperature properties can be deduced. By taking into account scattering of 
the quasiparticles the low temperature behaviour of the transport coefficients can 
be calculated, and equations for the coherent excitation of quasiparticle-quasihole 
pairs lead to a description of collective modes such as zero sound. Application and 
further discussion of the Landau phenomenological approach can be found in ref- 
erences [19] [20]. This approach was verified within the framework of many-body 
perturbation theory in the early 60s and detailed derivations of the basic theory can 
be found in the papers of Luttinger [6] and the books of Nozieres [7], and Abrikosov, 
Gorkov and Dzyaloshinskii [8]. Due to the mathematical complexity of the diagram- 
matic perturbation theory, however, some of the more intuitive ideas of Landau are 
lost. Here we intend to show that the renormalization group approach can make a 
useful conceptual link between the two approaches. For those not familiar with the 
renormalization group we very brief review of the philosophy of this approach as 
developed by Wilson for tackling the problems of critical phenomena, and also for 
the Kondo problem [17]. 

The renormalization group is a mapping R of a Hamiltonian H(K) } which is 
specified by a set of interaction parameters or couplings K = (Ki , K 2} . . .) into 
another Hamiltonian of the same form with a new set of coupling parameters K' = 
(K[ } K'21 • • •)• This is expressed formally by 

R{H(K)}=H(K'), (3) 

or equivalently, 

R(K) = K', (4) 

where the transformation is in general non-linear. In applications to critical phe- 
nomena the new Hamiltonian is obtained by removing short range fluctuations to 
generate an effective Hamiltonian valid over larger length scales. In the cases we 
consider here the transformations are generated by eliminating higher energy excited 
states to give a new effective Hamiltonian for the reduced energy scale for the lower 
lying states. The transformation is characterized by a parameter, say a, which spec- 
ifies the ratio of the new length or energy scale to the old one such that a sequence of 
transformations, generates a sequence of points or, where a is a continuous variable, 
a trajectory in the parameter space K. The transformation is constructed so that it 
satisfies 

R a ,{R a (K)} = R a+a ,(K). (5) 

The key concept of the renormalization group is that of a fixed point, a point K* 
which is invariant under the transformation, 

R a (K*) = K*. (6) 
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The trajectories generated by the repeated application of the renormalization group 
tend to be drawn towards, or expelled from, the fixed points. The behaviour of 
the trajectories near a fixed point can usually be determined by linearizing the 
transformation in the neighbourhood of the fixed point. If in the neighbourhood of 
a particular fixed point K = K* + <5K then, expanding i?„(K) in powers of <5K, 

R a {K* + 6K) = K* + L*JK + 0{6K 2 ), (7) 

where L* is a linear transformation. If the eigenvectors and eigenvalues of L* are 
O* (a) and A* (a), and if these are complete then they can be used as a basis for a 
representation of the vector <5K, 

SK = J2$K n O:(a), (8) 

n 

where 8K n are the components. How the trajectories move in the region of a par- 
ticular fixed point depends on the eigenvalues A*, if we act m times on a point K 
in the neighbourhood of a fixed point K* then, from (8), we find 

iC(K* + 6K) = K* + Y: SK n X: m O* n (a), (9) 

n 

provided all the points generated by the transformation are in the vicinity of the fixed 
point so that the linear approximation (9) remains valid. Eigenvalues with A* > 1 
are termed relevant and the corresponding components of 8K n in (9) increase with 
m while those with eigenvalues A* < 1, termed irrelevant, get smaller with m. Those 
with A* = 1 to linear order do not vary with m and are marginal. The eigenvalues 
of the linearized equation lead to a classification of the fixed points: stable fixed 
points have only irrelevant eigenvalues so <5K — > 0, unstable fixed points have one 
or more relevant eigenvalues and the trajectories are eventually driven away from 
the fixed point. A marginal fixed point has no relevant eigenvalues and at least one 
marginal one. Whether a trajectory is ultimately driven towards or away from the 
fixed point in this case depends on the non-linear corrections. 

The application of the renormalization group to critical phenomena is based 
on the fact that at a second order phase transition the correlation length becomes 
infinite so that at the critical temperature the system appears the same on all 
length scales. As a result at the critical point the parameter changes on applying the 
renormalization group transformation at the critical point are small and tend to zero 
as the length scale becomes large compared to the lattice spacing. This behaviour 
corresponds to the asymptotic approach of the trajectory to a fixed point of the 
renormalization group transformation. For small deviations away from the critical 
temperature, this scale invariance no longer holds and under the renormalization 
group transformation the system trajectory is eventually driven away from the fixed 
point. For this to occur the fixed point must be an unstable one. For critical 
behaviour of the form \T — T\ v the critical exponent v can be calculated from the 
relevant eigenvalues of the linearized renormalization group equations about this 
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fixed point (for further details see for example [21]). The universality of the critical 
exponents follows when a whole class of systems scale to the same fixed point, and 
consequently all have the same critical behaviour. 

In the situations we consider here the renormalization group transformations 
are associated with the elimination of higher energy states to generate an effective 
Hamiltonian for the lower lying levels. If the system does not undergo a phase 
transition, or develop an energy gap, then we expect the physics on the lowest 
energy scales not to change significantly so that we find a limiting form for the 
effective Hamiltonian as the energy scale is reduced to zero. This Hamiltonian 
must correspond to a stable or marginal fixed point of the renormalization group 
transformation. It will be this fixed point Hamiltonian, and the leading correction 
terms, which determines the lowest lying excitations of the system and hence the 
thermodynamic behaviour as T — > 0. 

As one universal form of low energy behaviour is that of a Fermi liquid it seems 
a natural question to ask whether we can reformulate Fermi liquid theory within the 
renormalization group framework. The idea of a Fermi liquid fixed point has been 
used frequently in condensed matter theory (see for example [22]) but very little work 
has been done in clarifying the concept and relating it to existing phenomenological 
and microscopic theory. Renormalization group calculations for fermion systems, 
in which higher energy excitations are eliminated and an effective Hamiltonian ob- 
tained for the lowest energy scales, are difficult to carry out. There have been a some 
calculations for translationally invariant systems and lattice models (see [23] [18] and 
references therein) but the most successful calculations of this type have been those 
for magnetic impurity models [17] [24], the Kondo and Anderson models for a mag- 
netic 3d transition or 4f rare earth ion in a metallic host. These are relatively 
simple systems yet with non-trivial physics, with a local moment regime at higher 
temperatures and a crossover to Fermi liquid behaviour at low temperatures and 
eventually to a spin compensated ground state. Apart from the extensive renormal- 
ization group calculations there are perturbational results [25] as well as an exact 
solutions for the thermodynamics [26] for these models. It will be instructive to look 
at the renormalization group results for these systems for which we have a rather 
complete quantitative description from a range of theoretical approaches. The in- 
sights gained from these specific examples will enable us to make conjectures, and 
devise methods, to apply to a very general class of systems. 

2 Fermi liquid theory of impurity models 

We base our discussion primarily on the Anderson model [27] as it gives a more 
general description of a magnetic impurity than the Kondo model because it includes 
the possibility of charge fluctuations at the impurity site. In its simplest form the 
model has an impurity d (f ) level e^, taken to be non-degenerate, which is hybridized 
with the host conduction electrons via a matrix element Vjj. When the interaction 
term U between the electrons in the local d state is included the Hamiltonian has 
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the form, 



k,<7 



k,<7 



where A(cu) = tt J2k |Fk| 2 <5(a; ~~ e k) is the function which controls the width of the 
virtual bound state resonance at in the non-interacting model ill = 0). In 
the limit of a wide conduction band with a flat density of states A(cu) becomes 
independent of u and can be taken as a constant A. For <C £f, td + U 3> ep, where 
ep is the Fermi level and | — ep|, | + U — ep| ^> A, the model is equivalent to the 
Kondo model with an antiferromagnetic interaction Jp = U A/Tr\ed — ep| \td-\-U — ep|, 
where p is the conduction electron density of states. 




Figure 1. A linear chain form of the Anderson impurity model with the impurity at 
one end (filled square) coupled by the hybridization to a tight-binding chain of conduction 
states (filled circles). Iterative diagonalization proceeds by diagonalizing a finite chain, 
adding a further conduction electron site to the chain, and repeating the process. 

We look first of all at the renormalization group calculations of Wilson who was 
the hrst to obtain a precise description of the low temperature behaviour of the model 
in the Kondo regime where the spin fluctuation effects dominate. We need not go 
into the technical details of the numerical renormalization group method he devised 
but it will help to have a brief overview to see what was involved. The model was 
cast in the form with the conduction electron states in the form of a tight-binding 
chain with the impurity at one end coupled via the hybridization (see figure 1). By 
iterative diagonalization of chains of increasing length, in which a finite number of 
the low lying states were retained at each stage, the lowest excitations on a decreasing 
energy scale were obtained (the details can be found in the original papers of Wilson 
[17], and Krishnamurthy, Wilkins and Wilson [24]). The transformation that maps 
the model on one energy scale to one on a reduced energy scale, when suitably 
scaled, is a renormalization group transformation. The higher energy excitations, 
which are eliminated at each stage, modify the couplings of the model. In general 
further interaction terms, not present in the original model, are introduced. It can 
be shown that the effective Hamiltonian for the lowest energy scales found by this 
technique can be expressed in the same form as the original Anderson model but 
with modified (renormalized) parameters, [28] [29], A = tt |Fk| 2 <5(a; ~~ e k)? U, and 
id, corresponding to the effective Hamiltonian, 

H eS = Y, e ^ C i,a C Xa + Y,(^la C k,a + K C l,a^ ( U ) 

k,<j k,<j ° 
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where the energy level ed, a is measured with respect to the Fermi level. 

For the local moment regime of the model there is only one energy scale involved 
which is the Kondo temperature Tk- The important point which makes the renor- 
malized effective Hamiltonian (f f ) tractable at low temperatures is that it describes 
excitations from the exact ground state so that the interaction term U between ex- 
citations only comes into play when there is more than one excitation from the 
ground state. From the renormalization group perspective the interaction terms 
are the leading irrelevant corrections to the free fermion fixed point. They tend to 
zero under the rescaling of the renormalization group as the low energy fixed point 
is approached (T — > 0). If a single particle excitation or single hole excitation is 
created this interaction term plays no role and can be omitted and the one electron 
Hamiltonian remaining can be diagonalized and written in the form, 

H eS (U = 0) = J2^clc^ (12) 

1,(7 

where e\°} are the one-electron energies, and c\ a and c t a the corresponding creation 
and annihilation operators. This is the fixed point Hamiltonian as T — > 0, and corre- 
sponds to free or non-interacting fermions. Comparing this with the Landau theory 
we can identify the excitation energy as the energy of the non-interacting quasi- 
particles in the free energy functional (I). Rather than diagonalize (II) explicitly 
for the excitation energies we can solve for the Green's function of the quasiparticles 
by the equation of motion technique which gives a closed set of equations for U = 0. 
This gives 

GWM = r 1 — x- (13) 

the quasiparticle energies correspond to the poles of this Green's function. The cor- 
responding spectral density or quasiparticle density of states pd{^) can be straight- 
forwardly be deduced and is given by 

7T [(to - e d ) 2 + A 2 ] 

corresponding to a Lorentzian resonance, the so-called Kondo resonance. 

For U ^ (11) can be written using the single particle eigenstates of (12) as a 
basis, 

#eff = e l^ C \,a C l,a + U a l '«*»«/"' ' C \,] C l> ,t4< ,l C l"> ,| ( 15 ) 

1,(7 1,1' 

where 

As it is the excitations of the interacting system from its ground state that 
are described by the effective Hamiltonian (11), it is appropriate to transform this 
Hamiltonian to operators which describe the single particle excitations, 

C L = P\,a C l,a = VI, a %a > , (17) 
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4,a = V C l,a = h \,a < £ F, (18) 

where the ground state is such that p/,<r|0) = and /i/ j(7 |0) = 0. The interaction 
term has to be normal ordered in terms of these operators so that H e ^\0) = 0, and 
the Hamiltonian describes interactions only between excitations from the ground 
state. In the mean field approximation the expectation value of operators in the 
interaction terms such as p] ^p t , ^p]„ ^p^,, ^ are approximated by 

(phPi;iPi»,lPi'",l) = (P^i^Pi^ipl'dPi"^ 6 ','' 6 '",''"- ( 19 ) 

The quasiparticle energy of the Landau theory in the presence of other excitations 
can be identified as the effective one particle energy in this approximation, 

eV = eg + U\m\ 2 J2 K| 2 (<W,-.)- (20) 

where (<*m// j(7 ) = {p]i a Pii a ) for q/ > ep and <*m// j(7 = —{h\i a hi, a ) for q/ < q?. as |ct/| 2 
is proportional to 1/N S , where N s is the number of sites, because the scattering 
potential is due to a single impurity, the energy shift in (18) is of the order 1/N S . 

As (p\ a Pi a ) — > as T — > the quasiparticle interaction does not contribute to 
the linear term in the specific heat which can be calculated from the non-interacting 
quasiparticle Hamiltonian (12). Hence the impurity specific heat coefficient 7i mp is 
given by 

7im P = M £ f)- (21) 

where pdiyj) is the impurity quasiparticle density of states defined in equation (14). 
in calculating the susceptibility the quasiparticle interaction has to be taken into 
account. Using the mean field approximation for this interaction given above the 
impurity susceptibility can be deduced and expressed in the form, 

Xi mP = ^y^^M(l + ^M) (22) 

for a flat wide conduction band. The total charge susceptibility at T = 0, \c = 
dn /dep } where n is the expectation value of the total number operator for the 
electrons, can be calculated following precisely the same argument. The result for 
the impurity contribution is 

Xc,im P = 2p d (e F )(l - Up d (e F )). (23) 

eliminating the term in U between (22) and (23), and the term in /^(ep) using (19), 
gives the well known Fermi liquid relation for the impurity model relating 7i mp , Xim P 
and Xc,imp? 

4ximp + Y . - 67imp (OA) 

(^b) 2 ^™" TT'fcg- 1 ] 

The interesting physics of this model occurs in the strong correlation regime when 
U is large and is well below the Fermi level so that the d electron at the impurity 
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is localized and has an occupation number ~ 1. This is the local moment regime 
in which the model maps into a Kondo model. At high temperatures T ^> Tk, 
where Tk is the Kondo temperature, there is a Curie law susceptibility associated 
with the local moment but with logarithmic corrections, ln(T/Tk), which reduce the 
effective moment. In the low temperature regime (11) is valid for T <C Tk when the 
impurity moment is fully screened giving a finite susceptibility at T — 0. The charge 
susceptibility on the other hand must tend to zero in the strong correlation limit as 
the d electron at the impurity is localized. The condition Xc,im P = is achieved in the 
quasiparticle picture by the interaction term U, which self-consistently constrains 
the many-body resonance at the Fermi-level so that the impurity occupation is 
maintained with = 1. This was hrst pointed out by Nozieres [30]. This condition 
can be used to deduce U . In the case of particle-hole symmetry 1^ = and the 
resonance in the quasiparticle density of states (14) is peaked at the Fermi level and 
Pdi^-p) = 1/ti"A. We then find on substituting into equations (22) and (23), 

U = ttA = 4T K , (25) 

where the Kondo temperature, which is the only relevant energy scale for the model 
in the local moment regime, is defined by Xim P = (5 , / u b) 2 /4Tk. 




0.0 1.0 2.0 3.0 4.0 

U/jiA 



Figure 2 A plot of the renormalized parameters U and A for the symmetric Anderson 
model in terms of the bare parameters U and A. In the comparison of these parameters 
with 4Tk for U 3* ttA the value of Tk is given by (27). 
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These parameters give for the 'x/7' or Wilson ratio i?, 



R 



4vr A; B Xim P Xim P /xo 



3(fl , /«B) 2 7imp 7imp/70 



(26) 



and so is enhanced over the free electron value of unity, where Xo and 70 are the 
susceptibilities and specific heat coefficients for the conduction electrons alone. The 
argument used here is essentially a reformulation of the one originally given by 
Nozieres [30]. The result R = 2 was first derived by Wilson from his numerical 
renormalization group results [17]. 

As the Anderson model is integrable and exact solutions exist for the thermody- 
namic behaviour it is possible to deduce the renormalized parameters exactly over 
the full parameter regime. For the symmetric Anderson model the parameters / ' 
and A have been calculated using the exact Bethe ansatz results [26] from equations 
(24) and (26). These are shown in figure 2 for the symmetric model (rid = 1) plotted 
as a function of U/ttA. At weak coupling U ~ U as one would expect and U and 
A are independent energy scales. At strong coupling there is only the single energy 
scale Tk given by 

T K = U (J^ V2 e -^/sA+,A/2U^ (2?) 



with U and A given by (25). Using these equations the form of the non-interacting 
quasiparticle density of states pd(u) given by (14), normalized by 1/ttA, for vari- 
ous values of U is shown in figure 3 illustrating the exponential narrowing of the 
resonance at the Fermi level in the Kondo regime (17 > 7rA) as U is increased. 



0.40 
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Figure 3. The quasi-particle density of states p~d(u) normalized by 1/trA as a function 
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of U for U/ir=0, 0.5, 1.0, 1.5 (in order of decreasing width) and A = 1.0 for the particle- 
hole symmetric Anderson model. 

This approach can be used for other magnetic impurity models, such as the N 
fold degenerate Anderson model (U = oo) and the n = 2S multi-channel Kondo 
model [28]. Though explicit renormalization group calculations do not exist for 
most of these models one can conjecture the form of the quasiparticle Hamiltonian 
to describe the Fermi liquid regime and deduce the Fermi liquid relations. These 
relations can be checked with the exact results known from either the Bethe ansatz 
or microscopic Fermi liquid theory. For these models there are more renormalized 
parameters to take into account but in the strong correlation regime they can all 
be expressed in terms of the Kondo temperature. We give an example in hgure 4 
where we plot the quasiparticle density of states in the Kondo regime for the A" fold 
degenerate model for different values of A" [28]. The case N = 2 corresponds to 
the Kondo resonance at the Fermi level of the large U Anderson model as shown in 
hgure 3. With increasing A" this resonance narrows and moves towards Tk above 
the Fermi level, becoming a delta function at u = Tk in the large A" limit. Much of 
the low energy and low temperature physics of this model can be deduced from this 
asymmetric form for the quasiparticle density of states. 

0.8 



0.6 



3 0.4 



0.2 



o.o 

-4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 




Figure 4. The quasiparticle density of states p(oj) for the N-fold degenerate Anderson 
model (U = oo) in the Kondo regime for N = 2,3,4,6,8, as indicated, as a function of 
cj/Tk where Tk is the Kondo temperature. 

What general conclusions about Fermi liquid theory might we be tempted to 
draw from these explicit renormalization group calculations? First of all, there 
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is strong parallel with the original phenomenological approach of Landau. The 
renormalization group calculations have given us an explicit Hamiltonian to describe 
the low energy excitations. The mean field expectation value of this Hamiltonian 
corresponds to the expectation value of the Landau total energy functional. The 
fact that this energy functional had to be expressed in terms of the deviations in 
the occupations of the states from their ground state values has its parallel in the 
renormalization group fixed point Hamiltonian, the Hamiltonian had to be normal 
ordered because it only describes the excitations from the ground state. The fact 
that the mean field theory is asymptotically exact as T — > is because the deviation 
in the occupation of these states tends to zero in this limit; this corresponds to a 
low density of excitations. For the mean field treatment of the Hamiltonian to be 
valid there must be no singular scattering between the quasiparticles at low energies. 
Such singular scattering occurs in one dimensional models, and in higher dimensions 
with attractive interactions. We will discuss this more fully later. 

A microscopic derivation of these Fermi liquid relations for the symmetric An- 
derson model has been given in a series of papers by Yamada and Yosida [25]. Their 
calculations were based on an expansion in powers of U for the symmetric model. 
Though it is not possible to sum explicitly the dominant terms in the strong correla- 
tion limit in this model exact relations can be deduced corresponding to Fermi liquid 
theory in terms of the local self-energy evaluated at the Fermi level {u = 0), 

its derivative and the irreducible four point vertex. Setting up a correspondence 
with the above we can identify the renormalized parameters, e<i, A and U, in terms 
of these [29]. With z, the wavefunction renormalization factor, given by 

1 (28) 



1 - E'(0)' 

where prime denotes a derivative with respect to u (evaluated at the Fermi level 
to = 0), the renormalized parameters are given by 

ed,* = z(ed,<, + M0,0)), A = zA, U = z 2 T tl {0, 0), (29) 

where r CT)CT /(u;,u/) is the irreducible four point vertex function. The self-energy has 
two arguments set to zero to indicate that it is evaluated not only for u = but 
also at zero temperature and in zero field. 

We have succeeded in showing that from the Wilson normalization group ap- 
proach for the impurity model a quasiparticle Hamiltonian can be derived corre- 
sponding to a renormalized Anderson model, and that this in the mean field ap- 
proximation gives the Landau total energy functional. Both the Landau and the 
renormalization group approaches give exact results as T — > 0. If we try to extend 
the renormalization group calculations beyond the Fermi liquid regime by gener- 
ating an effective Hamiltonian over a higher energy range then further and more 
complicated interactions will come into play as the renormalization group trajecto- 
ries are no longer under the exclusive influence of the low energy fixed point. On 
higher energy scales the effective Hamiltonian approach only becomes useful when 
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the trajectories once again become dominated by a particular fixed point. For the 
crossover regime the effective Hamiltonian approach is not feasible and we have to 
rely on explicit numerical renormalization group calculations of the energy levels 
for predictions. The spectral density Pd{^>) of the d electron Green's function over 
the full energy range for the symmetric Anderson model, as calculated by the nu- 
merical renormalization group [31], is shown in figure 5. The quasiparticle density 
of states only describes the small region u near the Fermi-level, the extra peaks at 
to ~ iU/2 are associated with the broadened 'atomic' peaks. The magnetic impu- 
rity problem is rather a special case in that we have a numerical renormalization 
results over all the relevant energy scales. As such calculations are difficult to carry 
out in general for other interacting fermion systems it is worth while to look at 
another renormalization procedure in which we make a perturbation expansion in 
terms of the fully dressed quasiparticles [32]. This corresponds to a reorganization 
of perturbation theory in close analogy with that developed originally for quantum 
electrodynamics where the expansion is in terms of the particles with their observed 
masses and charges (interaction strengths), ie. fully dressed. In the held theory case 
the reorganization of the perturbation theory was a necessary one to eliminate the 
ultraviolet divergences and obtain finite results. Due to the finite upper cut-offs in 
condensed matter physics problems such a reorganization of perturbation theory is 
not a necessary on one in order to get finite results, nevertheless, it makes a lot of 
sense to work with the fully dressed quasiparticles, particularly for systems where 
the renormalization effects are very large such as in heavy fermion systems. What 
is more it integrates the microscopic perturbation theory derivation of Fermi liquid 
theory with the more intuitive Landau approach. 



0.5 



o.o 



U/kA-- 


= 


U/7lA = 


= 1.5 


U/7lA = 


= 2 


U/7tA = 


= 4 


U/7lA = 


= 6 







-20.0 -10.0 0.0 10.0 20.0 

m/A 



Figure 5. The spectral density irApdiuj) for the d electron Green's function of the 
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symmetric Anderson model for values of U/ttA as indicated (see reference [31]). 



3 A renormalized perturbation approach 

We continue with the Anderson model and consider the retarded one particle double- 
time Green's function for the d-electron expressed in the form, 

Gda{t0) = \l - A f \ fTT' ( 30 ) 

u — 6^ — A(uj) + zA (u) — Zj a {uj) 
where A(cu) was defined earlier and A(cu) = |Vk| 2 /(iu — ejj). We continue to 

work in the wide band limit where A(cu) is independent of u and A(cu) — > 0. The 
function Y, a (u) is the proper self-energy within a perturbation expansion in powers 
of the local interaction U . We will need the corresponding irreducible four point 
vertex function Y a ^ a i (cu, u/), which is a special case of the more general irreducible 
four point vertex function r <7j(7 / j(7 // j(7 ///(a;, u/; u/', u'") with a" = cr } a 1 " = cr' } u" = u 
and lo'" = io 1 . Our aim is to reorganise this perturbation expansion into a more 
convenient form to consider the strong correlation regime which corresponds to a 
model with U large, at low temperatures, and with weak magnetic fields. 
Our first step is to write the self-energy in the form, 

E» = £,(0) + cuE^O) + £r», (31) 

which is simply a definition for the remainder self-energy E^ em (cu). Using this ex- 
pression the Green's function given in equation (2) in the wide band limit can be 
written in the form, 

G*,H = Z ~ f ■ (32) 

to - e d + zA - h a (uj) 

We assume that the general theorem of Luttinger [6], that the imaginary part of 
E(cu) vanishes as to 2 at u = 0, so that z is a real quantity. The 'renormalized' 
quantities, e, A, are defined by (29), and the renormalized self-energy E(cu) by 

t a (u) = zY?r(u). (33) 

The next step is to introduce rescaled creation and annihilation operators for the 
d-electron via 

C d,a = ^Z~ C \,a C d,a = ^~ C d,a- (34) 

We can now rewrite the Hamiltonian (1) in the form, 

H = Hqp — H c , (35) 
where will be referred to as the quasiparticle Hamiltonian which can be written 
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as -^qp^ • The Hamiltonian describes non-interacting particles and is given 

by 

17 k,<j k,<j 

and t[W is the interaction term, 

= Un dA n d , h (37) 

with U given by (29). The second term in equation (35) will be known as the counter 
term and takes the form, 

H c = XiJ2 ~ c \/ c d,o + X 2 h dA h d ^ (38) 
where Ai and A 2 are given by 

X 1 = z£(0, 0), A 2 = ^ 2 (r Tii (0, 0) - U) (39) 

Equations (36)— (37) are simply a rewriting of the original Hamiltonian (10), so 

equation (35) is an identity. 

We note that by construction the renormalized self-energy Y, a (u) is such that 

E„(0,0) = 0, E;(0,0) = 0, (40) 

so that Yj a (uj) = 0((jJ 2 ) for small u, on the assumption that it is analytic at u = 0. 
As r <7j(7 (0,0) = we also have 

rv(o,o) = u(i -8^,). (4i) 

To develop a theory appropriate for the low temperature regime we follow the 
renormalization procedure as used in quantum field theory so that we can make a 
perturbation expansion in terms of our fully dressed quasiparticles (see for instance 
reference [33]). We take our renormalized parameters e d} A and U as known and 
reorganise the perturbation expansion in powers of the renormalized coupling U . 
The full interaction Hamiltonian is — H c . The terms Ai, A 2 and z are formally 
expressed as series in powers of U, 

oo oo oo 

X 1 = J2 A^CT, X 2 = J2 ^2 ] U n , z = J2 z(n) U n - (42) 

n=0 n=0 n=0 

The coefficients A^, A^ and are determined by the requirement that condi- 
tions (40) and (41) are satisfied to each order in the expansion. The perturbation 
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expansion is about the free quasiparticle Hamiltonian given in equation (36) so that 
the non-interacting propagator in the expansion for the thermal Green's function is 

Gfffe) = , . r . v (43) 

%uj n - e d + zAsgn(cu n ) 

where u n = (2n + l)vr//3 and f3 = 1/ksT. The retarded function (30) is deduced by 

analytic continuation. 

We propose to make a direct contact with both the Landau phenomenological 
and the microscopic formulations of Fermi liquid theory. The occupation number of 
the d-level can be deduced from the Friedel sum rule [34] which takes the form, 
at T = and has the form 

ndA H ) = 2~k [ A J ' ^ ' 

in terms of the self-energy E(cu) (in a finite magnetic field H) for the expansion 

in powers of U for the 'bare' Hamiltonian (10). It takes the same form for the 
renormalized expansion, 

ndA H ) = 2~k [ A J ' ' 

This follows directly from (44) by writing it in terms of ed, CT , A and E from equation 

(33) as the common factor of z in the argument of (45) cancels. The quasiparticle 

interaction plays no role as T — > and H — > as E(0,0) = and in this limit 
nd,<r corresponds to the non-interacting quasiparticle number. As the effects of the 
quasiparticle interactions go to zero as T — > due to the cancellation with the 
counter term giving E(0,0) = and E'(0,0) = 0, the specific heat coefficient of 
the impurity 7i mp is due to the non-interacting quasiparticles so we obtain the free 
electron result (21). 

Other thermodynamic results for the low temperature regime can be obtained 
from the lowest order term in the renormalized perturbation for E, the tadpole 
diagram shown in figure 6(a), which gives 

£W(u, H } T) = U(n { %(^ H } T) - n$(0, 0, 0)). (46) 

There is no wavefunction renormalization to this order so = and also to this 

order f ^(cu, to 1 ) = U, \^ = and = Un^l(0,0). The complete cancellation by 
the counter term only occurs for H = and T = 0. The impurity spin susceptibility 
at T = of the impurity to first order in U can be calculated from gfj^ (rea,T — n d,|) /2, 
by substituting the self-energy from (46) into equation (45), and then differentiat- 
ing with respect to H } to give the earlier result (22), and similarly for the charge 
susceptibility (23). It is not obvious that the results to first order in U for these 
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quantities should be exact. However, there are renormalized Ward identities [32] [25] 
which can be derived from charge and spin conservation, 



dh 



dp 



(47) 



where h = gpsH/2. The spin and charge susceptibilities can be derived from the 
exact relations, 



Ximp 



(5^b) 2 



p d (0)(l-dt/dh) 



p d (0)(l + Up d (0)). 



(48) 



an 



d 



Xchnp = 2/5 d (0)(l + d^/dp) = 2^(0)(1 - Up d (0)). (49) 
using (45), (46) and (47), confirming that the earlier results. These relations show 
that there are no higher order contributions in the renormalized expansion in U for 
these quantities beyond hrst order. 




(a) 



(b) 



Figure 6. The two lowest order diagrams in the renormalized perturbation theory. 

Exact results for the impurity Green's function (32) to order u 2 follow from the 
calculation of the second order diagram for E shown in figure 6(b). There is no 
second order counter term contribution to this diagram from A 2 as = 0. There 
is a contribution to z which is required to eliminate the contribution from the linear 
term in u to this order and this is given by 



,(2) 



7T 



12) 



4tt 2 A 2 



(50) 



Calculation of this diagram to order to 2 gives 



Im S CT (u;, 0) 



U 2 -- 2 



2A(7rA) 



+ oV). 



(51) 



The spectral density of the non-interacting renormalized d-Green's function de- 
scribes the Kondo resonance and the u 2 terms in E CT (u;,0) give the low frequency 
corrections to this picture. 
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There is a temperature dependent contribution to E CT (0,T) to first order in U 
given by (46). For the particle-hole symmetric model, and more generally in the 
Kondo regime n d — > 1, this vanishes as nd j(J (0 } T) = ra<i )CT (0,0) = 1 and the leading 
order temperature dependence arises from the second order diagram 1(b). This gives 
the result, 

ImS ff (0, T) = ^ U J k B T2 + (T 4 ). (52) 
V ' ; 2A(vrA) 2 V ; ^ ' 

These results can be used to calculate the leading order T 2 contribution to the 

impurity conductivity cr- lin p(T) which is given by 

l + 2(i?-l) 2 ) + 0(T 4 )}. (53) 




A 

where R is the Wilson ratio given by 

R=l + U/ttA. (54) 

This conductivity result was first derived by Nozieres for the Kondo regime [30] and 
for the more general case by Yamada [25]. The formula for the Wilson ratio is a 
more general one than that given earlier and goes over to (26) R — > 2 in the Kondo 
regime on using the relation (25) between the renormalized parameters. Hence all 
the basic Fermi liquid results can be obtained within the renormalized expansion 
up to second order in U . The wavefunction renormalization factor z is required to 
relate the bare spectral density Pd{^>) to the quasiparticle density of states pdiyj) 
but does not enter explicitly the calculation of the thermodynamics which can be 
expressed entirely in terms of the renormalized parameters q, A and U . 

We see that the effects of the counter terms play no really significant role in 
the Fermi liquid regime. The renormalized perturbation theory with the counter 
terms, however, goes beyond the Fermi liquid regime. In this perturbation theory 
nothing has been omitted so that in principle calculations can be performed at high 
temperatures and high fields allowing the bare particles to be seen. If the primary 
interest is in these regimes then the more appropriate starting point is the model 
in terms of the bare parameters rather than the renormalized ones. It should be 
possible to estimate the renormalized parameters from the bare ones by a separate 
variational calculation for the single and two particle excitations from the ground 
state. In principle it is possible from the results of the renormalized perturbation 
theory to estimate the bare parameters from the renormalized ones by inverting (39) 
and (42), 

e d = (e d -X 1 )/z, A = A/z } U = (U-X 2 )/z 2 } (55) 

where Ai, A 2 and z are implicit functions of q, A and U . In the very weak coupling 

limit U <C 7rA z — > 1 and e d — > Q, A — > A and U — > U } as we can seen in 
the results for the symmetric model in figure 2. The most useful regime for the 
renormalized approach, however, should be in calculating corrections to Fermi liquid 
theory. The calculation with counter terms gives a systematic procedure for taking 
such corrections into account. 
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4 Renormalized Perturbation Theory for Trans- 
lationally Invariant Systems 

The renormalized perturbation approach described in the previous section gives us 
a very clear picture of the physics of the low temperature regime for a magnetic 
impurity as described by the Anderson model. It combines elements of the renor- 
malization group, microscopic and phenomenological Fermi liquid approaches, as 
well retaining many of the more intuitive elements of the original Landau theory. In 
this section we investigate whether it is possible to generalize the approach to a sys- 
tem which has translational invariance, and deal with the complications which arise 
when the self-energy depends on the momentum vector k as well as the frequency 
U. We consider a quite general model of fermions in states classified by momenta k 
and spin a described by a Hamiltonian, 

h = E e k,.4, ff c k, ff + \ E ^(q)4 1 +q, ffl C k 2 -q, ff2 C k 2 ^ 2 C k 1 ^ 1 ' ( 56 ) 

k,(T k-^ ,k 2 ,<7 2 

where = k 2 /2m with m the mass, and an interaction V(c\) which is a function of 
the momentum transfer q. It will be useful to express the Hamiltonian in a more 
general form, 

H = J2 e ^ c la c ^ + r Sl 2 ^3-4( k i' k 2 5 k 3,k4)4 3jCT3 4 4 ^ 4 c k2iCT2 c kiiCTi , (57) 

k,(T ks,(TS 

where the interaction is antisymmetric under the exchanges (ki,<7i) <->■ (k 2 ,cr 2 ) and 
(k3,cr 3 ) <->■ (k4,cr 4 ). For translational invariance and spin conservation we have 

ki + k 2 = k 3 + k 4 <Ti + cr 2 = cr 3 + cr 4 . (58) 

For a spin independent interaction depending only on the momentum transfer q as 
in(56) with V(q) = V(— q), then the antisymmetrized form is 

r Sl 2 ^ 4 ( k i 5 k 2,k3,k 4 ) = 6 auas 6 a2:ai V(k 3 - ki) - S aua J a2:as V(k 4 - k 4 ). (59) 

We can then write (56) as 

k,<7 

\ E r l% 2 ^ 3 -4( k l' k 2 5 k l + C l 5 k 2-q) C k 1 +q, ff 3 C k 2 -q, ff4 C k 2 , ff2 Ck 1 ,. 1 , (60) 
k-^ ,k2,q,(js 

The interacting one electron Green's function can be expressed in the usual form, 

GA^co) = 1 (61) 

to — e k — L(k, to) 
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where E(k,cu) is the corresponding self-energy. To follow our previous prescription 
we expand the self-energy in a Taylor's series to first order about the Fermi surface, 
and keep the remainder term, 

= E(Jfep, e F ) + (u - £F ) 9E( ^ £F) + (Jfe - k ¥ ) d ^ k e¥) + E-(k,c) (62) 

where we have assumed = e(|&|). We can rewrite this as 

G„(Ku) = . * (63) 

- e fc - h(k,u) 

with z, the wavefunction renormalization factor on the Fermi surface, given by 

1 

and with the renormalized energies given by 

h = z [{e k - e F ) + (k- k F ) mk d F k eF) + E(£ F , e F )) . (65) 

Similarly to (33) we define a renormalized self-energy by 

t{k,u) = zE rem (£,co), (66) 
with to measured relative to the Fermi level. We introduce rescaled operators, 

(67) 

and define a quasiparticle Green's function, 



G„(k,u) = , 1 ~ ( , v (68) 

to — £ k — E(fc, cuj 

The full irreducible renormalized four vertex is defined by 

f (7l(72(7 3 (74 (ki,cui,k2,cu2,ki + q,cui + u,k 2 - q,u 2 - u) = 

z 2 T rJirJ2rJsrJi (k 1 ,uj 1 ,k2,LJ2,k 1 + q,u 1 +u,k 2 -q,oj 2 -u). (69) 

So far every step has been an obvious generalization of the impurity case. When it 
comes to defining the renormalized interaction we have to be a little careful. Due 
to repeated particle-hole scattering the vertex has a singularity such that there is 
no unique way of taking the limits u — > and q — > 0. We leave this problem for 
the moment and assume that we can find a suitable antisymmetrized renormalized 
quasiparticle interaction (ki, k 2 , ki + q, k 2 — q) and hence construct a quasi- 

particle Hamiltonian, 

k,<7 

\ E KZ^^iMM + q,k 2 - q)4 i+q ^ 2 - q ,^ 2 W ( 7 °) 
k 1 ,k 2 ,q,(js 
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As earlier, the ground state is defined as the vacuum state with no excitations and 
the Hamiltonian has to be normal ordered and expressed in terms of quasiparticles 
and quasiholes so that the interaction only comes into play when more than one 
excitation is created from the (interacting) ground state. Following the steps we 
used in the previous section we write our original Hamiltonian (56) in the form, 

H = Hqp — H c , (71) 

where H c is the Hamiltonian for the counter terms and takes the form, 

#c = £(a<°> + a( 6 )(*-M)4,A,.+ 

k,<7 

7 E A lil 2 ^( k l' k 2' k l +q, k 2-q)Ck 1 +q,.3 £ k 2 -q,. 4 Ck 2 ,. 2 Ck 1 ,, 1 - (72) 



4, , 

ki,k 2 ,q,CTs 



where 



A(°) = / S( ^' £f) , \^ = z(Z(k F ,e F )-e F ) (73) 



an 



d 



Ali 2 ^( k i> k 2,ki+q,k 2 -q) = 

te(k ls k 2 , kx + q, k 2 - q) - * 2 r£> (k ls k 2 , k x + q, k 2 - q) (74) 



We can define a renormalized perturbation expansion, as earlier, in powers of V with 
the coefficients of the counter terms and the z factor chosen to satisfy the conditions, 

E(£ F , £F ) = 0, mk d F k 6F) =0, E'(fcF,e F ) = 0, (75) 

for zero field and at zero temperature. These conditions reflect the fact that the 
quasiparticle energies and Fermi level used in the quasiparticle Hamiltonian are al- 
ready fully renormalized. We now consider how to choose V^^ 2 (ki, k 2 , ki + q, k 2 — q) 
and how to fix the counter term (ki, k 2 , ki + q, k 2 — q). In the renormalized 

perturbation approach the renormalized interaction V^^ 2 (ki, k 2 , ki + q, k 2 — q) is 
a matter of choice, because it is added and subtracted, but there should be an op- 
timum choice for calculations of the low temperature behaviour as there was in the 
impurity case. An obvious choice is 

v;^(k 1} k 2 , k t + q , k 2 - q ) = f^ 2 ^(k 1} o, k 2 , o, k t + q , o, k 2 - q , o) (76) 

which is the quasiparticle scattering function. Its antisymmetry follows from its 
definition. When we come to evaluate the lowest order Hartree diagram of the 
renormalized theory, this has to be evaluated in the zero momentum exchange limit 
q — > 0. For a more general retarded u dependent interaction this diagram has to 
be evaluated in the limit q — > followed by u — > (there is a small imaginary 
contribution to u associated with the time ordering for this diagram). Repeated 
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quasiparticle-quasihole scattering gives a singular contribution to the retarded in- 
teraction and it is important to take this into account before we take the limits 
so that we can take them in the correct order. Similarly the Fock term, which is 
evaluated for q — > k 2 — ki followed by u — > 0, also has a singular contribution from 
quasiparticle-quasihole scattering. 



+ 



k+q 






k 1+ q k 2 k 1+ q k 2 k f q 



k, k f q 




k [+ q 



Figure 7. Particle-hole diagrams corresponding to the integral equation (77). Fhe 
unfilled square vertex corresponds to U^^(ki, k 2 , ki + q, k 2 — q) and the square filled 
vertex to f P"] 1 ^^ (ki, k 2 , ki + q,k 2 - q,u). 

If we denote the renormalized vertex corresponding to the sum the quasiparticle- 
quasihole series shown in figure 7 by F^"^ (ki, k 2 , ki + q, k 2 — q,cu) then we obtain 
the integral equation, 

f^U^ki, k 2 , kx + q, k 2 - q,cu) = U^(k 1} k 2 , kx + q, k 2 - q) 



+ E / UZ"{kiXM + q,k" - q)F(k",q;cu)f^^ 4 (k" - q,k 2 ,k",k 2 - q 



to 



J2 I ^"(ki, k", k 2 + q, k" - k 2 + kx - q)F(k", k 2 - kx - q; 



to 



f^* (k" ~ k 2 + ki + q,k 2 ,kx + q,k»Jk" 



where 



[/; 3 i ; 4 2 (kx,k 2 ,ki + q,k 2 -q) = 

^(kx, k 2 , k! + q, k 2 - q) - A^ 2ff3ff4 (k 1 , k 2 , k! + q, k 2 - q) 

and the propagator for the quasiparticle-quasihole pair F(k" } q;iu) is given by 

1 /(ek») - /(ek»-q) 



F(k",q;cu) 



F (27r) 3 iid + £ k » - £k"-q 



(78) 



(79) 
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where Vo is the volume. If we take the limit u — > in (77), then as (79) does not 
vanish in this limit it is clear that the counter terms must be non-zero as we have 
implicitly already taken some of these contributions to Y aia2a3ai (kx, k 2 , ki + q, k 2 — q) 
into account. We have to choose the counter term so that at T = we satisfy the 
condition expressed in equation (76). This implies 

V^(k 1} k 2 , kx + q, k 2 - q) = ^(ki, k 2 , kx + q, k 2 - q) 

+ E / 0"( k i> k "> k i + q> k " - q)^( k "> q; o)%(k" - q , k 2 , k", k 2 - q )jk" 
- E / £C'"( k i> k "' k 2 - q> k " - k 2 + k i + q)^( k ", k 2 - k i - q; o) 

a' a" 

y; 3 > (k" + kx - k 2 + q, k 2 , kx + q, k")dk" (80) 

For a given VJ 1 ^ (kx, k 2 , kx + q, k 2 — q) let us denote the solution of (80) for 
Vim (kx , k 2 , kx + q, k 2 - q) at T = by (k x , k 2 , k x + q, k 2 - q) . The antisymme- 

try of Ul^ll (kx, k 2 , kx + q, k 2 — q) with respect to the transformation (kx + q, cr 3 ) <->■ 
(k 2 — q, cr 4 ) (or equivalently (kx,o"x) <-> (k 2 ,cr 2 )) follows from equation (77). 

We can now show that in the limit q — >■ [/^^(kx, k 2 , kx + q, k 2 — q) can be 
identified as the Landau quasiparticle interaction function /^^(kx, k 2 ), which is the 
more general form than that introduced in (1) which takes account of the possibility 
of spin dependent interactions. 

The contribution to the direct term in the integral equation (77) has a pole 
arising from F(k" } q;iu) for small q and u. At T = and as q — > we have, on 
continuing to real frequencies (ico — > u), 

F(k",q;cu) = 6(\k"\ - k F ) n "'? VF (81) 

(uj — n" ■ qvp) 

where vp = kp/m and n" = — k"/|k"|. 

We can follow the standard treatment given in reference [8], chapter 4, and 
represent f ^ 2<73<7i (k x , k 2 , k x + q, k 2 - q, to) by x*i* 3 ( k i> q, cu)x^ 2(74 (k 2 , near the 
pole. If we keep only the dominant terms for small q and u we can cancel the 
common factor x^ 2(74 (k 2 , q, cu) and then write 

^1*3 ( n i) = 7 ^Xa 1 a 3 (k 1 ,q,uj), (82) 

[to - nx • qv F ) 

where nx is a unit vector in the direction of kx, so that for q — > we find 

(u - nx • qupK^m) = nx • q-^ J (k 1? k", k 1? k")v^ (n")iO", (83) 

where the integration is over the solid angle 0". Equation (83) corresponds to the 
Landau equation for the zero sound collective modes (see[8]) so it enables us to make 
the identification, 

/; 3 ^(kx,k 2 ) = ^; 4 2 (kx,k 2 ,kx,k 2 ). (84) 
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Figure 8. The Hartree-Fock diagram in the renormalized perturbation theory. The 
vertex corresponds to U^lf (ki, k 2 , ki + q, k 2 — q) and the double line represents the 
dressed propagator. There is also a counter term arising from the first term in (72). 

In the renormalized perturbation expansion we should include the more general 
vertex with the repeated quasiparticle-quasihole scattering given by (77) in evaluat- 
ing the Hartree (-Fock) diagram, figure 8, so that we take account of all the singular 
contributions in the u — > limit. We then get the correct results when we take the 
limits q — > (q — > k 2 — ki) followed by u — > 0. Evaluation of the Hartree-Fock 
diagram is equivalent to taking the mean field approximation on the quasiparticle 
Hamiltonian, 

k,<7 

\ E ^4 2 ( k l' k 2, k l + ^ k 2-q)c kl +q,.3 £ k 2 -q,. 4 Ck 2 ,. 2 C kl ,, 1 - 
k 1 ,k 2 ,q,(js 

(85) 

This gives the total excitation energy in the form, 

E tot = E gs + E e V<^k,, + \ E UZl^iMMMSn^Sh^ (86) 
k,<r ki ,k 2 ,(Ti ,(t 2 

when we have spin independent interactions. This equation can be identified with 
the Fandau equation (1) as (k l5 k 2 ) = (k l5 k 2 , k l5 k 2 ). 

We can choose to work with the quasiparticle Hamiltonian (85) with the interac- 
tion ~Ul\l 2 (ki, k 2 , ki + q, k 2 — q). This has the advantage that the mean field theory 
corresponds to the Fandau total energy functional, as in the impurity case, and it 
has no explicit interaction counter term when the low energy collective quasiparticle- 
quasihole excitations are taken into account. It satisfies the conditions for the opti- 
mum choice of the quasiparticle Hamiltonian. 

The two scattering functions, (k 1} k 2 , k t + q, k 2 - q) and U^(k 1} k 2 , k t + 

q, k 2 — q), are related by the integral equation (80) for T = with Vlll\ (ki, k 2 , ki + 
q, k 2 — q) — >■ Villi (ki, k 2 , ki + q, k 2 — q). It follows from the definition (76) that in 
the q — > limit VJ 1 ^ (ki, k 2 + q, ki,k 2 — q) is the forward scattering function for 
two quasiparticles in the Fandau theory which we denote by A^^(ki,k 2 ), 

A™(k ls k 2 ) = V;^(k 1} k 2 MM)- (87) 

If we take the limit q — > in equation (80), however, the equation does not reduce to 
the usual one (see [8]) relating A^^(ki,k 2 ) and fl\ll (ki, k 2 ), due to the presence 
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of the exchange term. In this respect this derivation differs from the standard 
derivations of the microscopic Landau theory. However, equation (80) would appear 
to be more satisfactory because it is applicable for general q and includes the singular 
contribution at q = k 2 — ki arising form the pole in the exchange term in (77) at 
to = n" • (k 2 — ki — q)vp. 

To summarize the situation so far: 

(i) . We identify the quasiparticle Hamiltonian with the interaction (ki, k 2 , ki + 
q, k 2 — q) as the Hamiltonian which determines the low energy excitations of a Fermi 
liquid. We conclude that the mean field approximation to this model corresponds to 
the Landau total energy functional. This Hamiltonian has non-forward scattering 
terms (q ^ 0), which do not enter the Landau theory, but these do not contribute 

at the mean field level. For low energy quasiparticle-quasihole pairs in the region of 
the Fermi surface there is only a rather restricted phase space for such scattering. In 
three dimensions with a spherical Fermi surface and states ki, k2 close to the Fermi 
surface being scattering into k3, k4, then k3 and k4 can only lie on or near the circle 
on the Fermi surface where the plane perpendicular to ki + k2, and which passes 
through ki and k2, intersects the Fermi sphere. When the scattering is restricted 
to states on the Fermi surface (ki, k 2 , k 3 , k 4 ) can be expressed as U^^(6 } (f>) 

where 6 is the angle between the wavevectors ki and k 2 and <f> the angle betwen 
the plane containing the vectors (ki,k 2 ) and that containing (k 3 ,k 4 ). There have 
been attempts to write down an effective Hamiltonian for interacting quasiparticles 
with q dependent scattering and parameters chosen so as to reproduce the Fermi 
liquid theory for forward scattering (see for instance [35]). However we can see from 
the theory developed here the limitations and inconsistencies of any such approach 
which does not take account of the full Hamiltonian and include the counter terms. 

(ii) . We conclude also that Hqp given by (85) is the Hamiltonian for the Fermi 
liquid fixed point, plus leading corrections, that we would obtain if we could carry 
out fully the sequence of renormalization group transformations to remove the higher 
energy excitations. This is in analogy with the impurity case where we could show 
this directly. The mean field approximation for this Hamiltonian, like that of the 
impurity fixed point Hamiltonian, gives the Landau Fermi liquid theory. The molec- 
ular field nature of the Landau phenomenological theory has been emphasized in the 
review of Leggett [19]; here we can give it a more precise interpretation within the 
framework of the renormalization group. 

This form for the fixed point Hamiltonian is in general agreement with the calcu- 
lations of Shankar [18] based on a perturbative renormalization of diagrams up to one 
loop level. The leading corrections to the free particle Hamiltonian for translation- 
ally invariant systems were found by Shankar to be marginal whereas the leading 
corrections in the Wilson impurity calculation, discussed in section (2), were the 
leading order irrelevant ones. This is just a reflection of the fact that for an im- 
purity problem we are interested in the impurities effects which, after summation, 
scale independently of the size of the system, and which therefore behave as 1/Vo 
where Vq is the volume of the system. 
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(iii). The leading order thermodynamics as T — > can be calculated either 
by normal ordering the quasiparticle Hamiltonian, ignoring the counter terms, and 
taking the mean field theory, or by self-consistently taking into account the first order 
Hartree (-Fock) diagram together with the counter term. As in the impurity case 
the counter terms play no role in the Fermi liquid regime other than subtract off the 
ground state expectation values. For calculations beyond the very low temperature 
or low frequency regime the full Hamiltonian (71) with the counter terms must 
be used. This way of reorganising the perturbation expansion closely parallels the 
renormalized prescription used in perturbative field theory. Here, the approach is 
rather more similar to that for the </> 4 field theory than quantum electrodynamics 
because we do not have the problems associated with gauge fields. In the field 
theory case there is a choice about which point to make the renormalization, and 
this freedom is exploited in setting up the Callen-Symanzik renormalization group 
equations. This freedom does not seem to extend to the situation considered here 
because, if we expand about any point other than one on the Fermi surface, the z 
factor will in general become complex, this would seem to preclude the setting up 
Callen-Symanzik type of renormalization group equations in this context. 

This reformulation of Fermi liquid theory as a renormalized perturbation the- 
ory links the renormalization group, microscopic Fermi liquid theory, and the phe- 
nomenological theory of Landau within a single theoretical framework. There should 
be scope for its further development and application to a range of physical systems. 

5 Breakdown of Fermi Liquid Theory 

Most metallic systems are not Fermi liquids at very low temperatures because at 
some critical temperature they undergo a phase transition to a magnetic, supercon- 
ducting, or some other form of ordered state. Nevertheless if this transition occurs 
at very low temperatures it can usually be analyzed as an instability within a Fermi 
liquid, such as in the BCS theory of superconductivity. In the case of supercon- 
ductivity the quasiparticle scattering becomes singular at the critical temperature 
T c signalling the on-set of the superconductivity if there is an attractive interaction 
between the quasiparticles in one of the particle-particle scattering channels, no 
matter how small. More commonly, as for most magnetic transitions, there is some 
critical value of the interaction strength between the quasiparticles, which has to 
be exceeded if the scattering between the quasiparticles becomes singular. If quasi- 
particles are assumed to exist in one dimension the scattering is singular whatever 
the magnitude (non-zero) or sign of the interactions and there is no (finite) critical 
temperature, indicating the breakdown of the quasiparticle concept in this case [36]. 
As mentioned earlier Haldane has shown that for one dimensional systems with re- 
pulsive interactions the behaviour at low temperatures can be described within the 
framework of a Luttinger model [14]. This would seem to preclude any description 
in terms of the renormalized perturbation theory that we introduced in the earlier 



27 



sections as we assumed the self-energy to be analytic at u = to give a finite wave- 
function renormalization factor. However, we show here that this problem can be 
circumvented, at least for the spinless Luttinger model, by making our expansion 
at finite u and combining it with a form of poor man's renormalization. The cor- 
rect form for the Green's function is obtained with the exponents of the singular 
terms evaluated perturbationally. A form of renormalized perturbation theory with 
counter terms can then be set up for the regular contributions to the self-energy. 
We give an outline of the approach. 

In the Luttinger model there are two branches of the conduction states, a and 
6, with Fermi points at kp and — kp, and energies relative to the Fermi energy of 
vp(k — kp) and vp( — k — kp) } respectively. The Hamiltonian of the model [10] is 



where Vij(q) is the interaction which is symmetric with respect to i and j. In this 
version of the model the number of electrons of each type are conserved. We take 
a cut-off of D on the momentum transfer q, \q\ < D and a cut-off on the bands 
of W so that \k — kp\ < W for the a band and \k + kp\ < W for the b band with 
W > D. The Luttinger model is an approximation for a one dimensional system 
of interacting fermions in which the dispersion of e(k) is linearized about the two 
points of the Fermi 'surface'. 

The corresponding Green's functions can be written in the form, 



where k a = k — kp and kj, = —k — kp. 

We now reduce the cut-off D by 8D and separate the self-energy into two parts, 
Y*i(uj,ki) = E 8 - rem ^(cu, ki) + Y*\ 5D \uj, ki), where Y*\ 5D \uj, ki) is the contribution from 
terms in which there is any scattering from the band edge region and E 8 - rem ^(cu, ki) 
is the remaining part. For the self-energy Y*\ 5D \uj, ki) we expand in powers of (u — 
ki) I D and keep the leading terms. This should be sufficient for calculating the 
behaviour in the region {u — ki) < D. If E| 5£ V, ki) = ai (8D) + bi(8D)(u - ki) + 
0(lo — ki) 2 then for this region we can write 



where Zi (8D) = 1/(1 - bi(8D)), %(u,ki) = ^•(<5 J D)Ef em) (cu, k) and ai (8D) is ab- 



sorbed a a renormalization of kp and the Fermi velocity vp. Initially we calculate 
E 8 - 51 ^ (cu, ki) to second order in perturbation theory and justify this later. We do not 
attempt a complete solution at this stage but just try to pick up the most important 
contributions. For the system with a reduced cut-off we can set up a new pertur- 
bation expansion for the renormalized self-energy E 8 '(cu, ki) using the renormalized 



H = ^u f (&-&f)c, 



k 







(90) 
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parameters. This can be continued as long as all the steps in the calculation are valid 
so generating a set of renormalization group equations. The first order perturbation 
terms only renormalize the chemical potential so we concentrate on the second order 
terms for the a branch, which corresponds to figure 6(b), taking Vi^(q) = so that 
the a propagator scatters only with a particle-hole pair in the b branch. This gives 

V (2)/ m 1 fdk' f vlb(q)dq f((k' b - q)) - f((k' b )) , ( , , f((u 

K j K, h o) = —\ — \ -. — — 7— — W?) + f{{k a - ?)), (91) 

Vy J Ztt J Z7T \ito n + Zq — K a ) 

where we have now scaled the ks to absorb the Fermi velocity vp } and n(q) is the 
Bose distribution function. We integrate over k' for T = to obtain 

V (2)/ j, n 1 [ 2 f w (% - kq)0(q) - 0(-q + k a )e(-q) 

K j K. M = — u / v a ,b{q)q d q r- — — n • 92 

{2ttv y Y J (iLo n + 2q - k a ) 

We take account of the momentum exchange terms that fall within the band edge. 
It helps at this stage if we distinguish the upper and lower band edges taking the 
running values to be D\ and — D 2 . The second order contribution due to scattering 
in the band edges is given by 

2 Usi), - *, 2 6b 2 - _ 72> - k a )6D 2 1 

I (u + k a ) + 2D 1 2D 2 -(u-k a )\ V ; 

where 7; = v\ 6 (D 8 ')/4(27tuf) 2 and we have specialized to the case k a > and con- 
tinued to real frequencies. To eliminate a step function we have also displaced D\ 
by the transformation to D\ — k a — > D\. We can set up scaling equations by ex- 
panding to order 6D/D } the most important one being that for the wavefunction 
renormalization factor z a which is given by 

D ~ dD t [D _ dD 2 



lnz a = - ,M D ^-L (94) 

■Ja\uj+ka\ L)\ Ja\w-k a \ L) 2 

where a is a suitable constant subject only to the condition a 1. Due to the 
different regimes of validity of the expansions used the lower limits differ in the two 
cases. 

There is also the possibility of a renormalization of the running effective inter- 
action 7; to consider. There are two second order corrections for the a part of the 
interaction vertex for the simple model with scattering between the two branches 
only; one arising from the diagram in figure 9 and the other from the z a factor which 
in the renormalized expansion is absorbed into the vertex, and similarly for the b part 
of the vertex. It is straightforward to show that these two corrections cancel, and 
so if we take v a ^(q) independent of q then 7; is a constant 7; = 7 = v\ 6 /4(27tuf) 2 . 
This result can also be derived from a Ward identity, based on the conservation of 
the two species of fermions by the Hamiltonian [37]. 
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Figure 9. The second order contribution to the renormalized a vertex for the Luttinger 
model with the momentum exchange q to the b particle-hole pair in the interval D — SD < 
\q\ < D . 

As 7 is constant we can integrate (94) which gives the result, 

za = a ^ + w,-t,|\ (95) 

The renormalized self-energy E^(cu, k a ) for remaining the range of q, values in the 
range — a\u — k a \ < q < a\u + k a \ } is calculated. This is not singular at u = k a 
because of the u dependent upper and lower cut-offs. The final result is 

G a ^k a ) = t+^p±zML , (96) 
D^{u - k a ) 

where D is independent of a for a ^> 1 to leading order in 7. The renormalizations 
of k a have been absorbed into a renormalized Fermi velocity. This result corresponds 
to the exact form for the spinless Futtinger model, with the exponent 7 given to 
leading order. There is no quasiparticle excitation at u = k a as z a goes to zero 
at this point. In this derivation we do not follow the usual scaling procedure of 
changing the cut-offs on all the momenta. We reduce the cut-off on the momentum 
transfer q only. In progressively reducing this cut-off we take into account the short 
range parts of the interaction at each renormalization step. 

The second order terms neglected in setting up the scaling equations can be 
taken in to account by writing the self-energy as the sum, S a (cu, k a ) = E( reg )(u;, k a ) + 
E^ smg )(a;, k a ), where is the sing ular E( sing )(u;, k a ) is the singular contribution we have 
just calculated and E^ reg )(a;, k a ) is the remaining regular part. Then we can apply 
standard perturbation theory for E a (cu,A; a ) and deduce E^ reg ))(a; n , k a ) exactly to 
second order from 

Ei reg )(cu, k a ) = E a (cu, k a ) - (ki - k a ) + (u - k a )(e-^ 2 -~ k °V D2 - 1) (97) 

The terms arising from Ti^ ins \u n , k a ) cancel the non-analytic terms in E a (cu, k a ) to 
second order v a ^ and so regularize the expansion to this order. There are differ- 
ent possible procedures for calculating higher order corrections depending on the 
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range of u of primary interest. If the main interest is in the region near u = k a 
then by introducing an effective coupling v a ^ by v\ b = 47(27TUp), and compensating 
counter terms, one can fix the exponent 7 and calculate the regular terms in powers 
of 7. The counter terms are then determined by the condition that they eliminate 
any further the singular terms in (97) order by order in 7. We can thus obtain 
the essential features for the Green's function for the spinless Luttinger model from 
this simple renormalization group approach. The Green's function for the Luttinger 
model can be obtained exactly from a complete summation of diagrams [39] [40] or 
via bosonization [41] [38]. Here, however, we are concerned with finding approximate 
techniques which can capture the essential physics of the model that can be gener- 
alized to tackle some of the interesting physical models which are not tractable by 
the exact approaches. An example in this context where this approach might be 
applicable is the model of two directly coupled Luttinger liquids [42] which has been 
put forward to explain the enhanced superconductivity of the high T c compounds 
[43] [44]. 

The two dimensional Hubbard model is the most basic model which has been 
put forward to describe the electrons in the CuO^ planes of the high T c materials. 
This model has the form, 

H= Y.kda^a + UY.^h (98) 

(ij)o- i 

where t 8 J - is the hopping matrix element for electrons in the hybridized orbitals built 
up from the copper 3d(x 2 — y 2 ) and the oxygen per state, with U as a short range 
repulsive interaction. Bethe ansatz [45] and conformal invariance methods [46] [47] 
provide a comprehensive picture for the one dimensional version of this model but 
so far there has been limited progress in understanding the behaviour of the two 
dimensional model despite much active work in this field. The low energy excitations 
of the one dimensional model can be analyzed in terms of the holons and spinons 
of the general Luttinger model (with spin). The low energy features of the one 
electron Green's function can be obtained from conformal invariance methods [46]. 
The non-doped high T c materials correspond to the insulating half filled Hubbard 
model. For an insulating state described by the Hubbard model there must be a 
Mott-Hubbard gap in the one electron spectrum. In one dimension a gap exists no 
matter how small the interaction U [45]. In models of higher dimensionality this 
is only likely to develop at some critical value of U (= U c ). The superconductivity 
arises with relatively small doping from this insulating state so it is of some interest 
to understand the transition from a Fermi liquid to a Mott-Hubbard insulator as 
U is increased to the critical value U c , and also to determine the nature of the low 
energy excitations away from half-filling for U > U c where the electrons are the most 
constrained. In this parameter regime, because of the limited phase space available 
for scattering, the correlation effects are likely to be at their strongest and non- Fermi 
liquid states might occur or transitions to magnetic order (at half filling and finite 
but large U the model is equivalent to an antiferromagnetic Heisenberg model). This 
difficult regime appears to be the most tractable in the infinite dimension version of 
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this model where it has been shown that the self-energy is purely local (independent 
of the wavevector k) and the model can be mapped into an effective impurity model 
with an additional self-consistency condition [48]. By formally integrating out all 
the states except those of a single site in a functional integral representation of the 
partition function the model can be mapped on to an impurity problem with an 
effective hybridization matrix element. The hybridization can be absorbed into the 
resonance width function A(cu) of the Anderson model (10). The Green's function 
for the non-interacting version of this effective Anderson model can be put in the 
form, 

= u - Cd - 1/tt / A(uj')/(uj - uj')duj' ^ 

which is equivalent to (30). The self-energy S(cu) of the impurity model can be 
formally expressed in the form, 

EH = F([/,GgH), (100) 

where F is a functional of G^\(uj) and corresponds to the complete sum of the 

diagrammatic perturbation series. This sum is completely specified, given G^\(to) 
and U, and the functional F is universal in the d — > oo limit. The local or on-site 
Green's function for the infinite dimensional model has the form, 

GM = E ^rr = / — ^W e > < 101 ) 

k u ~ e k — M^J ■' u ~ 6 ~ M^) 

where p c ( e ) is the density of states of the non-interacting model (for a tight-binding 
hypercubic model this has a Gaussian form). The self-consistency condition arises 
from the fact that this Green's function corresponds to that of the equivalent impu- 
rity, and hence has to satisfy the condition, 

s( w ) = G i - 1 H-(^H)- 1 . (102) 

The resulting equations are those for the Anderson impurity model but, instead of 
the resonance width function A(cu) being specified as an initial condition, it has 
to be determined self-consistently. As one is now dealing with an impurity model 
there is a good chance that a reasonably accurate solution can be obtained. The 
difficult step is to calculate the self-energy given U for any given form of G^\(uj). 
Many methods have been devised for calculating the self-energy and the spectral 
density for the Anderson model, such a perturbation theory [49], Monte Carlo [50], 
the non-crossing approximation [51], and the renormalization group [52]. These all 
have their advantages and disadvantages. The easiest to carry out in practice is the 
second order perturbation theory which has been found to work very well for the 
standard impurity model with particle-hole symmetry up to values of U into the 
strong correlation regime (where there are three well defined peaks in the spectral 
density corresponding to that shown in figure 5). Second order perturbational results 
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for the local density of states p(co) of the infinite dimensional model on a Bethe lattice 
are shown in figure f for a value of U greater than then conduction band width but 
less than the critical value U c . If this is compared with the corresponding density 
of states for an impurity hybridized with a wide structureless conduction band as 
shown in figure 5 one sees the tendency of a gap to open up and the central peak 
associated with the quasiparticle density of states to become isolated. The local 
density of states for the non-interacting version of this model has the single peaked 
semi-elliptical form. In the results for the real part of the self-energy shown in figure 
II for the same model one can see a large negative gradient developing at the Fermi 
level resulting in a small value for the wavefunction renormalization factor z. As 
U — >■ U c the self-energy becomes singular at u = so z — > and the quasiparticle 
peak of the Fermi liquid disappears. This basic picture of the breakdown of the 
Fermi liquid state is confirmed in the Monte Carlo, non-crossing approximation and 
renormalization group approaches. There are some complications arising from the 
fact that there is a another solution corresponding to an insulating state which 
persists below U c but this solution appears to be a unstable one for T = [53]. 
With several groups actively working on this problem we might expect to get soon 
a detailed definitive theory of the breakdown of Fermi liquid theory at the Mott 
transition for this model. 
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Figure 10. The local density of states for the infinite dimensional Hubbard model at 
half-filling for a Bethe lattice (total bandwidth 2D = 2.0 for the non-interacting model) 
with U = 2.7 as calculated from second order perturbation theory for the impurity model 
with the self-consistency condition (102). 
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Figure 11. The real part of the self-energy S(w) for the infinite dimensional Hubbard 
model with the same parameters as given in figure 10. The dashed curve corresponds to 
the self-energy in the atomic limit. 

Breakdowns of Fermi liquid have also been proposed for some strong correlation 
models away from half- filling and with dimensionality greater than one [4] [54] . In 
this regime the breakdown of Fermi liquid theory is more unusual as the resulting 
state can be a metal but lacking a Fermi surface in the usual sense as the imaginary 
part of the self-energy does not vanish, as in the recent theory of Herz and Edwards 
[4]. There is an impurity model which shows non- Fermi liquid behaviour, the n- 
channel Kondo impurity model where an impurity spin S is coupled equally to the 
spins of conduction electrons in each of the n channels. The thermodynamics of 
this model can be calculated exactly by the Bethe ansatz [55] and non- Fermi liquid 
behaviour is found in the 'overcompensated' regime n > 28. The low frequency 
behaviour of the local Green's function can be determined from conformal field 
theory [57]. This model has been put forward to explain the anomalous behaviour 
of some uranium heavy fermion materials [56] . The relevance of these various non- 
Fermi liquid theories to strong correlation systems, and in particular to the high 
temperature superconductors, is likely to remain a subject of lively debate in the 
field for sometime. The non- Fermi liquid theories of strong correlation systems is a 
suitable topic for a future survey. 
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